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Abstract 

Central issues of the Dirac constraint formalism are discussed in relation to the algorithmic methods of commuta- 
tive algebra based on the Grobner basis techniques. For a wide class of finite dimensional polynomial degenerate 
Lagrangian systems, we describe an algorithmic scheme of computation of the complete set of constraints, their 
separation into subsets of first and second class constraints as well as the construction of a generator of local sym- 
metry transformations. The proposed scheme is exemplified by considering the so-called light-cone Yang-Mills 
mechanics with an SU (2) gauge structure group. 
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1 Introduction 

Lagrangians used for the description of fundamental particles, such as electrons and pho- 
tons, as well as quarks and gluons have a degenerate Hessian functions. This rather un- 
usual property, compared to standard Lagrangian mechanical models, profoundly modifies 
the whole mathematical description of classical evolution. It demands the physical inter- 
pretation of constrained variables (e.g. longitudinal components of the electromagnetic 
potential) and also requires the generalisation of the canonical quantisation scheme. From 
the mathematical point of view the new element of the Hamiltonian description of a degen- 
erate Lagrangian system is the involution analysis of the differential equations of motion. 
Its pivotal ingredients in the generalized Hamiltonian dynamics pQ- [S] are realised in the 
form of the Dirac scheme to determine constraints. This is related to [HI [7J E] the formal 
theory of differential equations [9]. The process of the determining all the integrability 
conditions that can not be derived using only the algebraic operations with the existing 
differential equations is just the "reproduction" of constraints in the Dirac formalism. 
Having a complete set of constraints we are able to identify the set of "truly" dynamical 
equations for this involutive system and therefore finally provide a deterministic classical 
evolution of the physical observables and perform the subsequent quantization. 

Effective completion to involution of systems of differential equations needed in field the- 
ories represents a very complicated challenge requiring sophisticated computer-algebraic 
methods [10]. Similarly the generalized Hamiltonian formalism also needs an efficient 
algorithmisation and implementation in a proper computer algebra software. 

In the present paper we apply the most universal algorithmic tool of commutative 
algebra, the Grobner bases as the main algorithmic ingredient of the generalized 
Hamiltonian dynamics for degenerate mechanical models with polynomial Lagrangians. 
In [12] it was already suggested to use the Grobner bases for the computation and separa- 
tion of constrains for such models. The underlying Dirac-Grobner algorithm is based on 
the facility of the Grobner bases method to manipulate with a polynomial in the phase 
variables modulo constraint manifold, and, in particular, to check whether the polyno- 
mial vanishes on the manifold. In the present note we propose some further algorithmic 
improvements and extensions aiming at the computational realization of the Hamiltonian 
reduction of degenerate mechanical system possessing local symmetries. 

It should be noticed that constructive ideas of the involution analysis of differential 
equations combined with those from the Grobner bases technique have culminated in the 
concept of involutive bases [13J as a special type of Grobner bases providing the efficient 
involutive algorithms [14J for construction of the involutive as well as the reduced Grobner 
bases. 

The plan of this paper is as follows. We start (Section 2) with a brief description of 
the main issues in the Dirac constraint formalism that should be put into an algorithmic 
form suitable for effective calculations. In Section 3 the ways to achieve this goal for 
finite-dimensional mechanical systems with polynomial Lagrangians are described. Then 
(Section 4) we consider the so-called light-cone SU(n) Yang-Mills mechanics as an in- 
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teresting example of constrained model for which the first algorithmic issue of the Dirac 
formalism, namely, construction of the primary constraints, can be performed for arbi- 
trary n. The remaining algorithmic issues of the Dirac formalism are illustrated in Section 
4 for this model specified [15] to the simplest nontrivial structure group SU{2). Finally, 
in Section 5 some conclusions are presented. 

2 The issues requiring algorithmisation 

Here we sketch briefly the basic notions and definitions from the Dirac constraint formal- 
ism for a finite dimensional degenerate Lagrangian system and make a list of the main 
procedures requiring an algorithmic reformulation. 

Consider an n-dimensional mechanical system whose configuration space is R n and 
the Lagrangian L(q, q) is defined on a tangent space as a function of the coordinates 
q := qx, q 2 , . . . , q n and velocities q := q x , q 2 , . . . , q n ■ 

The Lagrangian system is called a regular one if the rank r := rank\\Hij\\ of the cor- 
responding Hessian function Hij := d 2 L/dqidqj is maximal (r = n). In this case the 
Euler-Lagrange equations 

d fdL\ dL 

aH - a" = 0> l<i<n (1) 



dt \dqiJ dq, 
rewritten explicitly as 

H " • - — 

11 13 dqjdqi qj dqi 

can be resolved with respect to the accelerations (q) and there is no hidden constraints. 
Otherwise, if r < n, the Euler-Lagrange equations are degenerate or singular. In this case 
not all differential equations ([T]) are of second order, namely there are n — r independent 
equations, Lagrangian constraints, containing only coordinates and velocities. Passing to 
the Hamiltonian description via a Legendre transformation 

Pi ■= 7T 2 
dqi 

the degeneracy of the Hessian manifests itself in the existence of n — r relations between 
coordinates and momenta, the primary constraints 

<P ( a ] (P,q) = 0, \<a<n-r. (3) 

Equations ([3]) define the so-called primary constraints subset Si. This definition is implicit 
and therefore the first algorithmisation topic is: 

Issue I: Find all primary constraints describing the subset Si . 

From ([3]) the dynamics is constrained by the set Si and by the Dirac prescription is 
governed by the total Hamiltonian 

H T := H c + U a <f>^ , (4) 
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which differs from the canonical Hamiltonian Hc(p, q) = ptqi — L by a linear combination 
of the primary constraints with the Lagrange multipliers U a . 

The next step is to analyze the dynamical requirement that classical trajectories remain 
in Ei during evolution, 

^ ] = {Ht,^}^0. (5) 
In the evolutional changes are generated by the canonical Poisson brackets with the 

total Hamiltonian (jlj) and the abbreviation = stands for a week equality, i.e., the right- 
hand side of ([5]) vanishes modulo the constraints. 

The consistency condition (jHJ), unless it is satisfied identically, may lead either to a 
contradiction or to a determination of the Lagrange multipliers U a or to new constraints. 
The former case indicates that the given Hamiltonian system is inconsistent. 

In the latter case when is not satisfied identically and is independent of the multi- 
pliers U a the left-hand side of (jSJ) defines the new constraints. Otherwise, if the left-hand 
side depends on some Lagrange multipliers U a the consistency condition determines these 
multipliers, and, therefore, the constraints set is not enlarged by new constraints. The 
subsequent iteration of this consistency check ends up with the complete set of constraints 
and/or determination of some/or all Lagrange multipliers. 

The number of Lagrange multipliers U a which can be found is determined by the so- 
called Poisson bracket matrix 

-M{<l> a , fa], (6) 

where E denotes the subset of a phase space defined by the all constraints including 
primary <fy£\ secondary <f>@\ ternary (f)%\ etc., constraints, $ := (<f>$ , (j$ , • • • , 4$ ) 

E : a (p,g)=O, l<a<k. (7) 

The co-rank s := k — rank(M.) of matrix M represent the number of first-class constraints 
•01 , ip2 j • • ■ 5 ips j linear combinations of constraints <p a 

V>a (P, q) =J2 C "/3 (P> l) ' ( 8 ) 

p 

whose Poisson brackets are weakly zero 

{MPiQ),MPi<i)} = 1<oc,(3<s. (9) 

The remaining functionally independent constraints form the subset of second-class con- 
straints . 

This method of constraints determination in the Dirac formalism represents the partic- 
ular form of completion of the initial Hamiltonian equations to involution] the generated 
constraints are nothing else than the integrability conditions [6J, [T], [8] . 

Therefore the second algorithmisation challenge can be formulated as 

Issue II: Determine all integrability conditions and perform their separation 
into first and second class conditions. 
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First-class constraints play a very special role in the Hamiltonian description: they 
provide the basis for generators of local symmetry transformations. The knowledge of a 
local symmetry transformation is important because according to physical requirement 
the physical observables are singlets under the gauge symmetry transformations. 

So the next important algorithmisation problem is 

Issue III: Construct the generator of local symmetry transformation and find 
the basis for singlet observables. 

The last problem has direct impact on the process of Hamiltonian reduction, that is a 
formulation of a new Hamiltonian system with a reduced number of degrees of freedom but 
equivalent to the initial degenerate one [2j [161 EZ] • The presence of s first-class constraints 
and r := k — s second-class constraints guarantees the possibility of local reformulation 
of the initial 2n dimensional Hamiltonian system as a 2n — 2s — r dimensional reduced 
Hamiltonian system (cf. [7j). 

Therefore, the final fourth algorithmisation challenge we formulate here as 

Issue IV: Construct an equivalent unconstrained Hamiltonian system on the 
reduced phase space. 

3 How the algorithm works 

Here we extend the main ideas of [TJ] and describe the algorithmic basics that can be 
used to solve the problems stated in the previous section. In doing so, we restrict our 
consideration to an arbitrary dynamical system with finitely many degrees of freedom 
whose Lagrangian is apolynomial in coordinates and velocities with rational (possibly 
parametric) coefficients^! L(q, q) G Q[<7, g]. 

3.1 Primary constraints 

For degenerate systems the primary constraints (J3J) are consequences of the polynomial 
relations (j2j). These relations generate the polynomial ideal in Q[q, q,p] 

I p , qA = Id(ur =1 fe - dL/dq,}) C Q[p, q, q] . (10) 
Thereby, primary constraints §3$) belong to the radical ^I p , q of the elimination ideal 

I p ,q = J P,g,g n Q[P> 0\ ■ 

Correspondingly, for an appropriate term ordering which eliminates q, a Grobner basis of 
I p>g (denotation: GB(I p q )) is given by [TT| [TS] 

GB(I M ) = GB(I MA )nQ\p,q}. 

2 Throughout this section we use some standard notions and definitions of commutative algebra (see, for example, |18|). 
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This means that construction of the Grobner basis for the ideal ( [101) with omitting ele- 
ments in the basis depending on velocities and then constructing of GB(^Jl Ptq ) allows us 

to compute the set of primary constraints. If GB(^Jl P:Q ) = then the dynamical system 
is regular. Otherwise, the algebraically independent set $i of primary constraints is the 
subset $i C GB(^J 'l p , q ) such that 

V0(p, q)e$ 1 : q) G* Id($ 1 \ {<f>{p, q)}) . (11) 

Verification of ( 1TTT) is algorithmically done by computing the following normal form: 
NF(<p,GB(ld(^t\{(P})). 

Therefore, all the computational steps described above admit full algorithmisation by 
means of Grobner bases. In addition, the canonical Hamiltonian H c (p, q) is computed as 
NF(p iqi -L,GB(I MA )). 



3.2 Complete set of constraints and their separation 

The dynamical consequences (j^J) of a primary constraint can also be algorithmically ana- 
lyzed by computing the normal form of the Poisson brackets of the primary constraint and 
the total Hamiltonian modulo GB(^Jl pq ). Here the Lagrange multipliers U a in are 
treated as time-dependent functions. If the non-vanishing normal form does not contain 
U a , then it is nothing else than the secondary constraint. In this case the set of primary 
constraints is enlarged by the secondary constraint obtained and the process is iterated. 
At the end either the complete set $ of constraints (171) is constructed or inconsistency of 
the dynamical system is detected. The detection holds when the intermediate Grobner 
basis, whose computation is a part of the iterative procedure, becomes {1}. 

To separate the set $ = {<pi, . . . , <pk} into of first and second class constraints the 
Poisson bracket k x k matrix M ([6]) is built. Its entries are computed as normal forms of 
the Poisson brackets of the constraints modulo a Grobner basis of the ideal generated by 
set $. 

To construct s :— k — r; where r = rank(Wl) first-class constraints as linear combina- 
tions (jSJ) of constraints d7J satisfying @ it suffices to find the basis P = {p±, . . . ,Pk~ r } 
of the null space (kernel) of the linear transformation defined by M. Every vector p G P 
generates the first-class constraint of form p Q <j) Q . 

Now consider the s x k matrix (pj) a composed of components of vectors in P and find 
a basis T :— {ti, . . . , t r } of the null space of the corresponding linear transformation. For 
every vector t G T the second-class constraint is constructed as t a cf) a . 
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Thus the constraints separation can be done using linear algebra operations with the 
matrix M alone. Together with the Grobner bases technique this implies full algorithmi- 
sation for computing the complete set of algebraically independent constraints and their 
separation ( Issues I and II of Section 2). 

3.3 Generator of local symmetry transformations 

The local symmetries are generated by first-class constraints (cf. jS]) but the presence 
of the second-class constraints makes realization of the symmetry transformations very 
subtle. To overcome some of these difficulties one can effectively eliminate the second 
class constraints by changing the initial Poisson bracket to the new Dirac bracket defined 
as 

{/, 9}d ■= {f, 9} - {/, Xa}C~p{x/3, 9} , 

where Xa (1 < ot < r) denotes the second-class constraints, and the invertible r x r matrix 
C Qj a is defined as 

C Q/ 3 := {Xa, Xp} ■ 

Since for an arbitrary function / it follows that {/, Xo}d — the second-class constraints 
can be set to zero either before or after evaluating a Dirac bracket. This last evaluation, 
modulo the constraint functions, can be performed algorithmically exploiting the Grobner 
bases. After elimination of all second-class constraints follow to the Dirac conjecture pQ 
the generator G of local transformations is expressed as a linear combination of all first- 
class constraints 

G=jte?$ ) + f (12) 

(3=1 7=/ci+l 

and its action on phase space coordinates (q, p) is given now with the aid of the Dirac 
bracket 

5qi = {G,qi} D , 5pi = {G,pi} D . 

In (fl"2l the coefficients Ep and are functions of time t and the first sum includes k\ 
primary first-class constraints while the second sum contains the all remaining first-class 
constraints. Not all of the functions e in (1121) are independent ones. Here we briefly state 
how following the method suggested in [19j one can extract the irreducible set of functions 
from the set of e. The total time derivative of the gauge-symmetry generator (1121) is given 
in terms of the Dirac bracket of G and the canonical Hamiltonian: 

f- § + (13) 

Since the set of first-class constraints is complete, the Dirac bracket in the right-hand side 
of ([131) is 

{(/>n,H c }v = PnAv (14) 
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The unctions p^ v can be algorithmically computed by using the Grobner bases method. 
To perform this computation one can use, for example, the extended Grobner basis 
algorithm [20]. Given a set of polynomials F = {fi, ■ ■ ■ , f m } C Q[p, q] generating the 
polynomial ideal ld(F), this algorithm yields the explicit representation 

9i = hijfj (15) 

of elements in the Grobner basis {gi . . . , g n } of this ideal in terms of the ideal generated 
by polynomials in F. Therefore, having computed a Grobner basis for the ideal gener- 
ated by the first-class constraints and the corresponding polynomial coefficients for the 
elements in the Grobner basis as given in (iT5|) . the coefficients p^ u are easily computed by 
reduction [TU [TBI 120] of the Dirac bracket in (TTTI) modulo the Grobner basis expressed in 
terms of the first-class constraints <fi u . Note that one can similarly compute the algebra 
of first-class constraints 

{(pa, 4>I3}d = Qaf3~f4>-y , 

if the structure functions ^ Q/ g 7 are polynomials in p, q. 

The generator of local transformation is conserved modulo the primary constraints 

f E = =► 0f + 4 1 + 4 2 Wf = • (16) 

Since, by their construction, the constraints <f>j® do not vanish on the primary-constraint 
manifold Si, the relations (fl6l) represent the following system of differential equations on 
the gauge functions and 

e® + ef Pfh + 4 2 V = , (h + 1 < 7 < s) , (17) 

where the index /3 runs from 1 to k\, 7 runs from k\ + 1 to s and the functions p^ u are 
projected on to the subset S x . 

Since the differential system (TTTj) is underdetermined, one can express the functions e^jp 
in terms of arbitrary functions s^\t) and their derivatives [19] . Since this last procedure 
is algorithmic, this completes the algorithmic construction of the generator of the local 
symmetry transformation. 

The above described algorithmic procedures have been implemented as a Maple package 
(currently for Maple 10), and this package was used to perform the computations presented 
in the next section. 

It is worth noting here that the remaining part of Issue III as well as Issue IV still 
require an algorithmisation. 

4 Light-cone Yang-Mills mechanics 

Now we discuss the application of the general scheme described above to a mechanical 
model originated from Yang-Mills gauge theory formulated on the light-cone under the 
assumption of spatial homogeneity of the gauge fields. 
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The standard action of Yang-Mills field theory with structure group SU(n) in four- 
dimensional Minkowski space M4, endowed with a metric 77 is 

S := -= ( txF A *F, (18) 

where go is a coupling constant and the SU(n) algebra valued curvature two-form 

F := dA + A A A 

is constructed from the connection one-form A. The connection and curvature, as Lie 
algebra valued quantities, are expressed in terms of the antihermitian algebra basis T a , 

A = A a T a , F = F a T a . a = 1,2, ...,n 2 - 1. 

The metric rj^s enters the action through the dual field strength tensor defined in accor- 
dance with the Hodge star operation 



*F fMU := \ J det(^) €^ ua p F af3 , with totally antisymmetric tensor e^ap. 

The light-cone version of the theory is formulated using the frame where two vectors 
e± := ^ (e ± e 3 ) tangent to the light-cone are combined with the orthogonal pair e k , k = 

1, 2 . The corresponding coordinates are usually called (see, e.g. [213) light-cone coordinates 
x n _ f x +^ x - ^ x ±^j 

x . — 1 — I x — I — x I . x . — x . /c — 1.2. 
y/2 V 7 



The non-zero components of the metric are 7/ + _ = r]_ + = —i]u = —^22 — 1 • The connec- 
tion one-form in the light-cone basis is given as 

A :— A + dx + + A- dx~ + A k dx k . (19) 

By definition, the Lagrangian of light-cone Yang-Mills mechanics follows from the cor- 
responding Lagrangian of Yang-Mills theory if one supposes that the components of the 
connection one-form A in fTTTJ]) only depend on the light-cone "time variable" x + 

A±=A±(x + ), A k = A k (x + ). 

Substitution of this ansatz into the classical action (1181) defines the Lagrangian of light- 
cone Yang-Mills mechanics 



L:=— 2 (F«_ F a + _ + 2 F a +k F\ - F« F&) , (20) 

where g is the "renormalized" coupling constant g 2 = g 2 / (Volume) and the light-cone 
components of the field-strength tensor are given by 

f)A a 

Fl_ := — — + f abc A b ,A c _, 
ox + 
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8A a 

77>a A; | rabc Ab ac 

+k '~ dx+ + k ' 

F\ := f abc A b _ A c k , 

Ft ] :=i ahc A\A), i,j,k = 1,2. 

Therefore, (|2"01 determines the SU(n) Yang-Mills light-cone mechanics as 4(n 2 — 1)- di- 
mensional system with configuration coordinates A± , A k evolving with respect to the 
light-cone time r := x + . 

The Legendre transformation 

n„ := — — = , 
dA% 

dL 



dA a _ g 
dL 1 



f7 V 



:= — =- = — i abc A b Al 



" ' dAi g 2 



gives the canonical Hamiltonian 

He = y t^tt" - i abc A\ (At v ~ + At vr a fc ) + ^F&F& . (21) 

The non-vanishing Poisson brackets between the fundamental canonical variables are 

{A a ± ,n±} = 5 a b , {Aiy b } = 5 l k 5l. 



The Hessian of the Lagrangian system (120]) is degenerate, det 1 1 gj^j^ \ \ = 0, and as a 
result there are primary constraints whose computation by the algorithm of Section 3.1 
gives 

^ -=< = 0, (22) 
Xt ■= g 2 K + f abc A b _A c k = . (23) 

The non-vanishing Poisson brackets between these constraints are 

{Xhx b j}= 2gH abc A c _5 ij . 

According to the Dirac prescription, the presence of primary constraints affects the dy- 
namics of the degenerate system. Now the generic evolution is governed by the total 
Hamiltonian 

H T :=H c + U a (T)^ + V k a (r)x a k , 

where the Lagrange multipliers U a {r) and V k {r) are unspecified functions of the light-cone 
time t. Using this Hamiltonian the dynamical self-consistence of the primary constraints 
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may be checked. From the requirement of conservation of the primary constraints <p)p it 
follows that 

= = {< , H T } = P bc (A b _K + A b k rc k c ) . (24) 
Therefore, there are three secondary constraints 

<p?> := f abc (A b _7r; + A b k ^) = (25) 
which obey the SU(n) algebra 

{^ 2) ,^ 2) }=Uc^ 2) . 

The same procedure for the primary constraints xt gives the following self-consistency 
conditions 

= x a k = {x a k ,H c }-2gH abc V b A1. 

A further issue, the identification of the first class constraints among the primary con- 
straints x%i depends on the rank of the structure group. Below we specify to the simplest 
special unitary group of rank one. 

4.1 The SU(2) structure group 

Here we present the results of our computations performed for the case of 577(2) algebra 
where the structure constants are given by the totally antisymmetric three dimensional 
Levi-Civita symbol, i abc = e abc . 

Constraints and their separation. Computation of the complete set of constraints, as 
described in Section 3.2, gives nine primary constraints (f£, Xk an d three secondary con- 
straints y 9 i 2 \ in accordance with (!22l) and ( !25l) . Performing the separation of the primary 
constraints (l23il we find two additional first-class constraints 

A :=A a _ X a k , 

and four second class constraints 

(Atxl) A - 



{Al) 2 + (AL) 2 + (A 3 ) 2 ' 

The new first class constraints tpi are abelian, {if)i,i[)j} = 0, and also have zero Pois- 
son brackets with all other constraints, while the second class constraints xt± have the 
following non-zero Poisson bracket relations 

{x^,X b ^} = 2g 2 e abc A15 tJ , 

W ( a ] ,X b k± } = e abc X c k± . 

Summarizing, there are 8 functionally independent first-class constraints (p^^k-,^^ 
and 4 second-class constraints xt±- 
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Generator of local symmetry transformations. The presence of two first class constraints 
ipi raises the question of the existence of new local symmetries as well as the expected 
SU(2) gauge symmetry. To clarify this point we construct the corresponding generator 
of local symmetry transformation following Section 3.2. We start from the expression 

G = E4 1 Vi 1) + E^ + E4 2 Vi 2) , (26) 

a=l i=l a=l 

with the eight light-cone time-dependent functions e^(r), £%>{t) and r)i(r) , then com- 
pute the functions p (see eq. (JT4l)). Equation ( Tl6l) reads now as 

(4 2) + 4 1) -^c4 2) ^-^-)^ 2) = °- 

Therefore expressing in terms of the functions , the generator of local transforma- 
tion takes the final form 

G = (-4 2 ) + e^Al + m Af) 0« + ^ + 4 2 )0i 2 ) . (27) 

Analyzing the changes of the canonical coordinates A a and n a generated by (1271) we find 
that the abelian subgroup of the 5-parameter local symmetry is in some sense "inherited" 
from the rigid conformal symmetry of initial Yang-Mills theory. But now, instead of 
the conformal symmetry, the light-cone SU(2) Yang-Mills mechanics has the SL(2, R) 
dynamical group of symmetry. Moreover, the group action is accompanied by the abelian 
transformations generated by two constraints ipi . A detailed discussion of this symmetry 
realization will be given elsewhere. 

Hamiltonian reduction to unconstrained system. Now that we have the generator of local 
transformation, we can address the question of finding a set of suitable coordinates part 
of which represent the invariants of these transformations. Solving this problem will let 
us project our system onto the constraint manifold and thus determine the unconstrained 
Hamiltonian system. We refer for details to [15] , and here present the set of corresponding 
singlet variables (as an example of the solution of the second part of (Issue III ). We also 
give a result of subsequent implementation of a Hamiltonian reduction (Issue IV) of the 
"redundant" degrees of freedom associated to the symmetries generated by constraints 
<P ( a \ Va ] and VV 

Let us pass to a matrix notation: the 3x3 matrix A a b whose entries of the first two 
columns are A® and the third column is composed by the elements A a _. Now one can verify 
that the elimination of local degrees of freedom associated with the three constraints ip^ 
can be achieved by using the polar representation [22\ 

A = OS 

where S is a positive definite 3x3 symmetric matrix and the orthogonal matrix O is 
parameterized by three Euler angles. It turns out that these three angles represent the 
pure gauge degrees of freedom corresponding to the constraints p^ . 
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To find the gauge degrees connected with the remaining two abelian constraints ipi , ip2 
one can pass to a principal axes representation for the symmetric matrix S 



with the orthogonal matrix -R(xi?X2,X3) given in terms of the Euler angles (xi,X2,X3)- 
Now again it turns out that the two angles Xi an d X2 are pure gauge degrees of freedom. 
Solving for the remaining second class constraints Xi± leads to an unconstrained system 
which represents a free particle or, considering the complex solutions to the second class 
constraints, to a more interesting model, the so-called conformal mechanics. In this case 
the diagonal variable qi and the angular variable X3 together with the corresponding 
conjugate momenta p\ and p Xs are two unconstrained canonical pairs and their dynamics 
is governed by the reduced Hamiltonian 



which is a projection of the canonical Hamiltonian (1211) to the constraints shell. Finally, 
noting that p X3 is a constant of motion, the Hamiltonian f[2"S]) coincides with the Hamil- 
tonian of conformal mechanics with the coupling constant P X3 /^ . 

5 Concluding comments 

In this paper we have raised several issues for a constrained mechanical systems which 
require computational realization. We described how using the Grobner basis technique 
the computation and separation of the complete set of constraints as well as the con- 
struction of the local gauge transformations can be achieved in degenerate mechanical 
models whose Lagrangians are polynomials in coordinates and velocities. The remaining 
challenges, namely, the construction of a basis for singlet (gauge-invariant) variables as 
well as the subsequent Hamiltonian reduction still needs algorithmisation. However, a 
first step in this direction also has been performed. In systems with first-class constraints 
the configuration space should be factorized by the local symmetry group in order to find 
a gauge invariant basis. The infinitesimal structure of a local symmetry group is encoded 
in the generator of gauge transformations, and we have shown that its construction allows 
an effective algorithmisation. 

As an example of the effectiveness of the proposed algorithms light-cone Yang-Mills 
mechanics with the SU{2) structure group was analysed in details: we determined and 
separated constraints, constructed a local invariance transformation and found the equiv- 
alent unconstrained Hamiltonian system. 

For the SU(2) light-cone mechanics the computations with our implementation in Maple 
10, which is an improved and extended version of that given in [12], takes about 1 minute 
on a machine with a 1.7 GHz processor. This uses the standard Grobner package in the 
Maple library. Unfortunately, recent extensions of the Maple Grobner bases facilities with 



S = R T diag (gi , q 2 , <? 3 ) R 




(28) 
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the packages Gb and Fgb developed by J.C. Faugere [23] do not improve on the stan- 
dard package. Gb is slower for our problems while Fgb cannot deal with the parametric 
coefficients. For the same reason we cannot use our software GINV [21] to implement 
the involutive algorithms [H] for involutive or/and Grobner bases. Manipulation with 
parametric coefficients is essential for the Dirac formalism due to the presence of physical 
parameters (e.g. masses, coupling constants) in the initial Lagrangian, the Lagrange mul- 
tipliers in the total Hamiltonian (jl]) and the time-dependent functions in the generator 
(112j) of local symmetry transformations. 

Consideration of light-cone mechanics with n > 3 is under current study. Here we note 
only that a recent paper [25] on geodesic motion on the SU (3) group provides us with a 
useful parametrization suitable for this investigation. 
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